Quantum decoherence reduction by increasing the thermal bath temperature 
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The well-known increase of the decoherence rate with the temperature, for a quantum system 
coupled to a linear thermal bath, holds no longer for a different bath dynamics. This is shown by 
means of a simple classical non-linear bath, as well as a quantum spin-boson model. The anomalous 
effect is due to the temperature dependence of the bath spectral profile. The decoherence reduction 
via the temperature increase can be relevant for the design of quantum computers. 
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The rapid development of nanotechnology has opened 
the possibility of assembling materials and devices at the 
atomic scale, such as nanomaterials and nano electro- 
mechanical systems (NEMS) 1]. Because of quantum 
mechanical and statistical effects, the properties of these 
systems can be considerably different from their macro- 
scopic counterpart. Electrons in nanofabricated quantum 
dots as well as nanomechanical devices are possible candi- 
dates for quantum bits at nanometer scale, as suggested 
by recent experiments [2J. Quantum computers could 
take advantage of quantum superpositions in order to 
perform highly efficient parallel calculations. The main 
difficulty in the implementation is represented by deco- 
herence, due to the coupling with the environment 
This effect increases with the number of quantum bits 
and is devastating even with a small thermal bath cou- 
pling. Thus, in order to reduce the decoherence of quan- 
tum superpositions, a sufficiently small thermal distur- 
bance and environment coupling are required. The ther- 
mal disturbance is expected to be reduced by decreas- 
ing the bath temperature. This is true when the quan- 
tum system is coupled with a high number of degrees of 
freedom and the bath can be described as a linear noise 
source. On the contrary, at the nanometer scale, the 
quantum system can be directly coupled to few modes, 
whose power spectrum does not necessarily grow with 
temperature. Thus, it is possible to have a reversed ef- 
fect, i.e., the reduction of decoherence by increasing the 
thermal bath temperature. 

A similar counterintuitive behavior occurs for exam- 
ple in thermal rectifiers [3], where the heat flux can be 
suppressed by increasing the temperature difference of 
the junctions. The rectifier consists of two or three mu- 
tually coupled layers, whose phononic bands depend on 
the temperature via non-linear effects. As the ends of the 
chain are coupled with two heat baths at different tem- 
peratures T\ and T2, the overlaps among the phononic 
bands in each layer become functions of T\ and T2 and 
are not identical when the thermal baths are exchanged. 
As a consequence, the thermal flux can change when the 
rectifier is reversed. 

In this paper, we show that the decoherence rate of 
a quantum system coupled to a non-linear thermal bath 
can decrease with the temperature. As in thermal recti- 



fiers, this effect is due to the non-linearity which makes 
the shape of the bath spectrum dependent on the temper- 
ature. For our purpose, we consider a classical non-linear 
bath model, consisting of a one-dimensional particle in a 
double-well potential, strongly coupled to a linear bath. 
Furthermore, we also consider a full quantum mechan- 
ical bath model, consisting of a boson bath coupled to 
a two-state system, which is the low temperature limit 
of a double- well system under suitable conditions [5j. A 
probe quantum system is coupled to the non-linear bath 
and its decoherence rate is evaluated. 

The equation of motion of a one-dimensional particle 
coupled to a white noise thermal bath is mx = —V'(x) — 
mrfi + v '2Drj(t) , where m, V(x), 7 and D are the particle 
mass, the external potential, the dissipation and diffusion 
coefficients, respectively. The function i](t) is a Gaus- 
sian noise with (r](t)j](t')) — 5(t — t'). The fluctuation- 
dissipation theorem establishes that D = jmkbT, fcf, and 
T being the Boltzmann constant and the thermal bath 
temperature, respectively. We consider the double-well 
potential V(x) = —ax 2 /2 + bx 4 /A. By means of the 
rescaling x -> (a/bf^x, t -> (m/a) 1 / 2 t and T ->■ £-T, 
the equation of motion becomes 



bk h 



= x - x 3 - -fix + y / 2"/ 1 Trj(t), 



(1) 



where 71 = 7(m/a) 1 / 2 . The independent parameters are 
the rescaled dissipation coefficient and temperature. 

At finite temperature, the particle has a finite prob- 
ability to jump from a well to the other one. The hop- 
ping rate is @ R = (-\/27r7i) -1 exp[— (471T)" 1 ] (Kramers 
rate). In the limit of very low temperature, the hop- 
ping rate becomes very small and the particle remains in 
one well for a large time. As a consequence, the spec- 
trum of x(t) has a high zero- frequency peak. As the 
temperature is increased, the hopping rate between the 
two wells grows, the central peak of the spectrum reduces 
and two smooth peaks around the frequencies ±7ri? are 
generated. This behavior is shown in Fig. [TJ where we 
report the spectrum I(w,T) = 2Re ' (x{0)x(t)) T e lu>t dt 7 
for 71 = 0.4 and some values of the temperature. I(uj, T) 
has been evaluated numerically over a time interval of 
200, in the non dimensional units, and 50000 ensemble 
realizations. We have verified that the time interval is 
sufficiently large to ensure a negligible time cut-off of 
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the integral. In Fig. [2l we plot 1(0, T) (solid line) and 
5r 



[(Q/2)&i + ea 3 A x(t)] \if>), 




FIG. 1: Power spectrum of x(t) for some values of T. The 
vertical lines are placed at ttR for each spectrum. 

1(0, T)/T = K(0, T) (dashed line) as functions of T. The 
decreasing behavior of K(0,T) suggests a simple model 
of thermal rectifier, as will be discussed in an extended 
paper. Here we will use the decreasing behavior of 1(0, T) 
to show the anomalous behavior of decoherence induced 
by the non-linearity. We consider a two-state quantum 
system (A) coupled to the double- well oscillator (B), as 
sketched in Fig. [3] The overall system i?-(linear bath) is 
the non-linear bath seen by A. The linear bath and the 
anharmonicity of B guarantee that the overall system is 
a reservoir and a non-linear system, respectively. 




FIG. 2: 1(0, T) and K(0, T) as functions of the temperature, 
in logarithmic scale. 




FIG. 3: Schematic representation of the non-linear thermal 
bath, coupled to the quantum system A. B is the anharmonic 
oscillator. 

This classical bath model is suitable for a sufficiently 
high thermal bath temperature T. A full quantum ap- 
proach will be introduced later on. We consider an in- 
teraction Hamiltonian linear in x and assume that the 
quantum state lip) = 4>-i\ ~ 1) + "01 1 1) satisfies the 



Schrodinger equation z^ijp = 

where a\ are the Pauli matrices in the orthonormal ba- 
sis {| — 1), |1)}. n is the energy difference between the 
bare energy levels [|1) ± | — l)]/y/2 and e is the coupling 
coefficient between the quantum and classical systems. 
The variable x satisfies the classical Eq. (JIJ. We have 
seen that the spectrum I(u>,T) of the anharmonic os- 
cillator is a decreasing function of temperature only for 
sufficiently small values of \u>\. If the quantum system 
has a dynamics sufficiently slower than the classical hop- 
ping rate (f2 <C R), then it is sensitive only to the low 
frequency components of the bath spectrum. We sat- 
isfy the condition Q <C R by assuming that the two lev- 
els are nearly degenerate. For the sake of simplicity we 
set Q = 0. If \ip(0)) = (| - 1) + \l))/V% is the quan- 
tum state at time t = 0, then at a subsequent time t 
the density matrix associated with a single realization of 

/ x e -2ief*x(t)dt \ 

x(t) is p[x] (f ) = 2 I e2je p o x(t)dt 1 I and it is 

a functional of the thermal bath variable. This depen- 
dence is removed by averaging over every realization of 
x(t). The reduced density matrix is < p[x] > x = Pa and 
the modulus of its off-diagonal elements 



| (e 8 * J? "W*) | = C(t). 



(2) 



is a measure of the coherence between the superposed 
states | — 1) and |1). C is a good indicator since | — 1) 
and |1) are two perfect pointer states of a pure dephasing 
process. It ranges between the values (no coherence) 
and 1 (complete coherence). 

We will show that the decoherence rate of A is pro- 
portional to 1(0, T), provided that e is sufficiently small. 
Let t c be the correlation time of x(t) and St ^ t c be 
a finite time interval. Let us consider the finite differ- 



ence 



2C(t + 5t)-C(t)=C(t) \(e 2l ^l 



v(t)dt\ 



1 



For a 



large temperature range, we can assume (x 2 ) ~ 1, thus, if 
2e5t -C 1, then we can approximate the previous equation 



as C(t + St)-C(t) ~ -2e 2 C(t) f t +St dt' J^° l dt(x(t')x(t)) . 
Since St 3> t c , we can replace the second integral with 
1(0, T) and have C = -2e 2 I(0,T)C = -D(T)C. The 
decoherence rate D(T), being proportional to 1(0, T), is 
a decreasing function of T (see Fig. [5]). This equation 
is correct when the integration step St is much smaller 
than l/(2e) and much larger than t c , thus the condition 
l/t c ^> 2e has to be satisfied. For 71 = 0.4, the correla- 
tion time is about 10 at T ~ 0.25 and decreases for higher 
temperatures. At T ~ 1-2, t c is for example about 5. 
Thus, e has to be smaller than about 0.1-0.05. This clas- 
sical model demonstrates that the quantum decoherence 
rate is a decreasing function of the temperature when 
the quantum system is resonant with bath spectral com- 
ponents depleted by the thermal energy supplying. The 
non-linearity is a fundamental ingredient, since for linear 
baths the spectrum grows linearly with temperature at 
any frequency. 



r t+st 
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Our classical model works when high order energy lev- 
els of the double- well system are activated, i.e., for a 
sufficiently high temperature T. Now, we consider the 
opposite limit and assume that only the two lowest en- 
ergy states are populated. This is possible if the energy 
difference A between the ground and first excited states 
is much smaller than the higher transition energies. A is 
the quantum tunneling frequency between the two wells. 
Let the Pauli matrix a B be the position operator, whose 
eigenstates | ± 1) correspond to the particle in the left or 
right well. This quantum system (B) is coupled to a bo- 
son bath (spin-boson model) and to the two-state probe 
system (A) described the Pauli matrices a\. The overall 
system B- (boson bath) is the non- linear thermal bath of 
A (Fig. El). 

More precisely, we study a model with the Hamiltonian 
H = %a\ + T a h + z&a°b + Efc ( & k +a{) + u k a\a k ] , 
where hk and a\. are the annihilation and creation oper- 
ators of the bath oscillators. In the present model x is 
replaced by the Pauli matrix a B . In order to understand 
how the decoherence of A depends on its bare frequency 
n, we have to evaluate the two-time correlation function 
(^bW^bW) m the Heisenberg representation for e = 0, 
as done with the first model. For rj sufficiently small 
(see below for a more precise condition), we can use the 
rotating-wave approximation and find the master equa- 
tion for B (section 10.5 of Ref. 0) 



dp B 
dt 



A 

-i~2 Wb,Pb] +C-PB, 



where 



Cp B = -76" [°B°'bPB + PB<7 B °B ~ 2 ^bPB^b] 

-j b {h + 1) [a+a B p B + Pbo-b^b ~ 2 ^bPb^b] 



(3) 



(4) 



and a B 



(a B ±ia B )/2. 7b is a constant coefficient which 
depends on rj and the density of the boson modes at the 
frequency uj k — A and n(T) = [exp(A/T) — 1] _1 is their 
expectation value (atflfe). This equation is well-known 
since it describes the decay of a two-level atom coupled 
to optical modes. 

From the master equation ([3|) and using the quantum 
regression theorem, we have that the correlation function 
at equilibrium is [7| 



A ^> 1/t c . It is clear that for u> ^ A, the spectrum ini- 
tially increases with temperature, reaches a maximum, 
then decreases and tends to zero for T — > oo. In the 
particular case lu = A, I(u>,T) is a monotonic decreas- 
ing function for any value of temperature. With this in 
mind, we choose the frequency fl of system A in order to 
maximize the effect of temperature on the decoherence 
reduction. The system A has to be resonant with the 
part of the spectrum T) which is strongly depleted 
by the temperature increase, thus we assume f2 = A. 
This case simplifies the analysis and is sufficient for our 
purposes. 

At this point, we take e ^ and write out the master 
equation for A and B. If e <C l/r c , then it is obtained 
from Eq. ([3]) by adding the bare Hamiltonian of A and the 
interaction Hamiltonian to\o B to the bare Hamiltonian 
of B. In the rotating-wave representation (interaction 
picture), the overall bare Hamiltonian (Aa\ + Aa B )/2 is 
removed and the following replacements are performed: 



A.B 



iAt 



A,B^ 



(6) 



The condition e <C l/r c justifies a rotating-wave approx- 
imation also for the interaction between A and B. Thus 
we obtain the master equation 



dp 
dt 



ie[a+a B +a A a B ~,p} + £p 



(7) 



In order to eliminate the parameter A in n, we define the 
rescaled temperature T = T/A. 
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FIG. 4: Numerically evaluated coherence C as a function of 
time for some values of temperature T. 



,(o)4(t)) 



-\t\/r c (T) 



-iAt„ T+e iAt n 1 



(5) 



and n\ = 



(7 n (Tn> = 



where n T = (a+a B ) = 2co e sh(A / 2T) — -n = \»b"bi 
1 — nj are the probability weights for the energy lev- 
els of B and r c (T) = [{2n + l)7b] _1 is the temperature- 
dependent correlation time. The corresponding spectral 
function is I(u,T) = 1+T ^ C "-A)" + i+ j^+A)^ ■ which 
has two peaks centered in ±A, whose areas and widths 
are constant and linear functions of temperature, re- 
spectively. The rotating-wave approximation is suitable 
when the two peaks have a negligible overlap, i.e., when 



In this model, the vectors |1) and | — 1) are not perfect 
pointer states because of the presence of the bare Hamil- 
tonian of the system A, thus we can not use the previous 
definition of coherence. We need a generalized measure 
of the deviation from the unitary evolution. A good in- 
dicator is provided by the quantity C — y/2Ti[p A ] — 1, 
Pa being the reduced density operator of the system A. 
For pure states it is equal to 1 and tends to zero for com- 
pletely incoherent high temperature thermal states. It is 
easy to prove that this definition is equivalent to Eq. © 
when |1) and | — 1) are equally populated pointer states. 
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At thermal equilibrium, we have C = tanh[l/(2T)]. Ob- 
viously, it is a decreasing function of temperature and for 
T — > 0, oo takes the maximal and minimal values 1 and 0, 
respectively. The transient behavior of C is numerically 
evaluated by setting the probe system A in the pure state 
[| — 1) + |l)]/v2 at the initial time t — 0, with the system 
B thermalized at the temperature T. In Fig. [4] we report 
C(t) as a function of time for 76 = 1, e = 5 • 10~ 2 and 
some values of T. Note that, for T = 80, the decoherence 
rate becomes very small with respect to the value at low 
temperature. 

The master equation Q has been obtained by means 
of the conditions e <C l/r c <C A. Under the constraint 
e <C l/r c , it is possible to find an analytical relation, since 
some components of p{t) can be adiabatically eliminated. 
The following approximate solution for the reduced den- 
sity operator of the probe system holds, 

PA = Pth + e-^ia&l + c 3 & A ) + e- 2e2T " t c 2 a 1 A , (8) 

where pth is the thermalized density operator, Ck are 
constant parameters which depend on the initial condi- 
tion and e 2 r c = e 2 7~ 1 tanh(l/2T) = is a decreas- 
ing function of temperature. Notice that this solution 
could be obtained directly by the correlation function |(SJ| . 
By means of the rotating-wave approximation and the 
standard hypotheses used to derive Eq. (J3J) and the re- 
sult of the previous model (second order expansion in 
the coupling strength and Markov approximation), one 
finds that pa is solution of a master equation identi- 
cal to Eq. but with 7/, replaced by e 2 r c tanhi = 

e 2 7~ 1 ^tanh . Equation {5J is its exact solution. 
This derivation will be discussed in an expanded paper. 

As in the classical model, we find that the decoherence 
and relaxation rates are decreasing functions of tempera- 
ture. In this case the effect is much stronger. For T> 1, 
the decoherence rate is inversely proportional to the tem- 
perature and for a sufficiently high thermal noise it can 
be practically reduced to zero. The other interesting fea- 
ture is the reduction of by increasing the coupling 
parameter 7&. This effects are due to the depletion of 
bath spectrum I(w,T) at the frequencies resonant with 
the quantum system. 

There is an interesting link between the anomalous 
behaviors in the last model and the Zeno effect. As- 
sume that the two systems A and B are in the state 
I — 1)a|1)b at time t. Because of their mutual cou- 
pling, at a subsequent time t + At, they evolve to the 
state I — 1)a|1) b + eAt\l)A \ — l)s and the probability of 
the transition | - 1)a\1)b -> \1)a\ - 1}b is (eAt) 2 . In 
the meantime, the coherence of this superposition is de- 
stroyed by the thermal bath, which acts as an observer of 



the pointer states | — 1)b and Thus, if r c is the cor- 

relation time of B, the probability of the transition after 
a time At > r c is 1)^(1 _i) ~ [et c ] 2 ^ = e 2 T c At 
and we obtain heuristically the relaxation rate in Eq. l[8"]). 
Alternatively, the anomalous effect can be explained as 
due to a rapid change of the phases of | — 1)a\1)b and 
— l)s which reduces the coherence of the two states 
and consequently their transition probability. It is inter- 
esting to note that the decoherence rate suppression for 
T — ► co in the last model can be obtained by a condition 
established in Ref. [f| [equation (25) in that paper]. By 
means of this, one finds that the Hilbert space of A is 
decoherence-free, provided that n >■ 1. 

In conclusion, we have consider two examples of non- 
linear baths weakly coupled to a quantum system A and 
shown that the decoherence rate of A is a monotonic de- 
creasing function of temperature. This anomalous behav- 
ior is due to the temperature dependence of the spectral 
profile in non-linear baths. It is important to realize that 
at a microscopic scale non-linear effects can be relevant, 
as it occurs with thermal rectifiers. For more realistic 
models, the decoherence rate could not be a monotonic 
decreasing function of temperature. However our anal- 
ysis indicates that the decoherence rate is not necessar- 
ily minimal at T — 0, as for linear baths. Decoherence 
reduction by the thermal energy increase can have rele- 
vant consequences for the design of quantum computers 
at nanometer scale. An experimental test of this effect 
is feasible with the present technology and can be ob- 
served in nanometer devices [H, Q. For example, A and 
B could be a superconducting qubit and a microwave 
resonator [§]. This work was supported by Ente Cassa 
di Risparmio di Firenze under the project "Dinamiche 
cerebrali caotiche" . 
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